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We study the enhancement of cooperativity in the atom-light interface near a nanophotonic waveg- 
uide for application to quantum nondemolition (QND) measurement of atomic spins. Here the 
cooperativity per atom is determined by the ratio between the measurement strength and the deco- 
herence rate. Counterintuitively, we find that by placing the atoms at an azimuthal position where 
the guided probe mode has the lowest intensity, we increase the cooperativity. This arises because 
the QND measurement strength depends on the interference between the probe and scattered light 
guided into an orthogonal polarization mode, while the decoherence rate depends on the local in- 
tensity of the probe. Thus, by proper choice of geometry, the ratio of good to bad scattering can 
be strongly enhanced for highly anisotropic modes. We apply this to study spin squeezing result- 
ing from QND measurement of spin projection noise via the Faraday effect in two nanophotonic 
geometries, a cylindrical nanofiber and a square waveguide. We find, that with about 2500 atoms 
using realistic experimental parameters, ~ 6.3 dB and ~ 13 dB of squeezing can be achieved on the 


nanofiber and square waveguide, respectively. 


I. INTRODUCTION 


Cooperativity is a measure of the entangling strength 
of the atom-light interface in quantum optics. Originally 
introduced in cavity quantum electrodynamics (QED), 
the cooperativity per atom, C1, can be expressed in terms 
of the ratio of the coherent coupling rate to decoherence 
rates, C1 = g°/(TeT A) where g is the vacuum Rabi fre- 
quency, ľe is the cavity decay rate, and T4 is atomic 
spontaneous emission rate out of the cavity [1]. Alter- 
natively, we can write C1 = (o9/A)F, where oo is the 
resonant photon scattering cross section of the atom, A 
is the cavity mode area, and F is the cavity finesse. Ex- 
pressed in this way, cooperativity is seen to arise due to 
scattering of photons preferentially into the cavity mode, 
compared to emission into free space, here enhanced by 
the finesses due to the Purcell effect. Strong coupling dy- 
namics seen in pioneering experiments in atomic cavity 
QED [2, 3] is now a mainstay in quantum information 
processing in systems ranging from quantum dots [4— 
6] to circuit QED [7, 8]. The N4 atom cooperativity, 
Cn = (Na4o0/A)F = (OD)F, where OD is the resonant 
optical depth. In this configuration, the collective de- 
grees of the atom can be manipulated by their common 
coupling to the cavity mode. 

Cooperativity also characterizes the atom-light inter- 
face in the absence of a cavity. In free space, an atom 
at the waist of a laser beam will scatter into the for- 
ward direction at a rate k x (a9/A)7., where ys is the 
photon scattering rate into 47 steradians [9]. Here the 
single atom cooperativity can be expressed to be pro- 
portional to ratio of these rates, C1 «x K/Ys X o0/A. 
The N4-atom cooperativity, in a plane wave approxi- 
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mation, ignoring effects of diffraction and cloud geome- 
try [9], Cn x Naoo/A = OD. To be self-consistent, here 
the beam area must be very large, so C1 is very small, e.g. 
Cı ~ 10~®, but for a sufficiently large ensemble, the OD 
can be large enough to lead to entanglement between the 
collective atomic degrees of freedom and the light. In 
that situation, measurement of the light leads to back 
action on the ensemble, and for an appropriate quantum 
nondemolition (QND) interaction, results in squeezing of 
the collective spin [10, 11]. QND measurement-induced 
spin squeezing have been observed in free space dipole- 
trapped ensembles [12-14] and in optical cavities [15- 
17]. The rate of decoherence is set by the rate of optical 
pumping, Yop X Ys, and we can characterize the cooper- 
ativity/atom by Cy = &/op. 

In recent years nanophotonic waveguides have emerged 
as a new geometry that complements cavity QED, and 
can lead to strong cooperativity [18-24]. Notably, the ef- 
fective beam area of a tightly guided mode can be much 
smaller than in free space and propagate for long dis- 
tances without diffraction. As such, ¢9/A can be orders 
of magnitude larger than in free space, e.g., 79/A ~ 0.1, 
and contribute collectively for a modest ensemble of a few 
thousand atoms trapped near the surface of the waveg- 
uide. Moreover, in some cases the Purcell effect can fur- 
ther enhance forward scattering into the guided mode 
when compared with scattering into free space. Taken 
together, these features make nanophotonic waveguides 
a promising platform for the quantum atom-light inter- 
face. 

In this paper we show that one can achieve an addi- 
tional enhancement to the cooperativity in a nanopho- 
tonic geometry that is not possible in free space. In par- 
ticular, we consider the QND measurement of the collec- 
tive spin of an atomic ensemble via a Faraday interaction 
followed by polarization spectroscopy. In this configu- 
ration the polarimeter effectively performs a homodyne 


measurement, where the probe is the “local oscillator” 
that interferes with the light scattered into the orthogo- 
nally polarized guided mode [9]. This signal thus depends 
on the spatial overlap of the two orthogonal polarization 
modes at the position of the atom. In contrast, decoher- 
ence due to photon scattering into unguided 47 steradi- 
ans occurs at arate ys is determined only by the intensity 
of the probe. The net result is that the cooperativity per 
atom, C1 X «/7s, primarily depends on the strength of 
the orthogonal polarization mode, and this factor can be 
enhanced, especially for highly anisotropic guided modes. 
Counterintuitively, we will see that the strongest cooper- 
ativity arises when the atom is placed at the position of 
minimum intensity of the azimuthally anisotropic probe 
mode where the intensity of the initially unoccupied or- 
thogonal mode is maximum. 

We study the enhanced cooperativity for two nanopho- 
tonic geometries: a cylindrical nanofiber formed by ta- 
pering a standard optical fiber with cold atoms trapped 
in the evanescent wave, as recently employed in a variety 
of experimental studies [18, 25-38], and a nanofabricated 
suspended square waveguide, currently investigated at 
Sandia National Laboratories [39]. For each geometry 
we study the use of the Faraday effect and polarimetry to 
perform a QND measurement of the magnetic spins [40], 
and thereby induce squeezing of collective spins of ce- 
sium atoms. A dispersive measurement of the number of 
atoms trapped near the surface of an optical nanofiber 
was first performed in [41], and quantum spin projection 
noise was recently detected using a QND measurement 
with a two-color probe in [31, 38]. Previously, we studied 
QND measurement-induced spin squeezing mediated by 
a birefringence interaction [24]. We will see here that, 
through the enhanced cooperativity, QND measurement 
via the Faraday effect can lead to substantial squeezing, 
greater than 10 dB in some geometries, for 2500 atoms. 

The remainder of the paper is organized as follows. In 
Sec. II, we lay out the theoretical description of the QND 
measurement and the relevant measurement strength. In 
addition, we describe how decoherence is included in the 
model through a first-principles stochastic master equa- 
tion, here for the case of alkali atoms, and cesium in par- 
ticular. From this we will see how cooperativity emerges 
as the key parameter that characterizes the squeezing. 
We calculate in Sec. III the squeezing dynamics for the 
different nanophotonic waveguides, atomic preparations, 
and measurement protocols. We conclude with a sum- 
mary and outlook for future work. 


II. QND MEASUREMENT AND 
COOPERATIVITY 


The theoretical framework describing the propaga- 
tion of light guided in a nanofiber and interacting with 
trapped atoms in the dispersive regime was detailed in 
our previous work [24]. We review the salient features 
here and include the generalization to the square waveg- 


uide. 

For waveguides that are symmetric under a 7/2 rota- 
tion around the z (propagation) axis, there are two de- 
generate polarizations for each guided mode and for each 
propagation direction. Assuming a nanophotonic waveg- 
uide that supports only the lowest order guided mode, 
and restricting our attention to modes propagating in the 
positive z-direction, we denote uy(r,) and uy(r_,) as 
the horizontally and vertically polarized modes that adi- 
abatically connect to x and y linearly polarized modes, 
respectively, as the cross section of the waveguide be- 
come large compared to optical wavelength. Note, in 
typical nanophotonic geometries these guided modes also 
have a nonnegligible z component. For a cylindrically 
symmetric nanofiber, these are the well-studied HE41 
modes; for a square waveguide, these are the quasi-TEo, 
and quasi-TMp; modes, shown in Fig. (1). One can 
solve for the guided modes of a cylindrical fiber ana- 
lytically [24, 42, 43]. We use a vector finite difference 
method to numerically solve for the guided eigenmodes 
of the square waveguide [44] with core material of Si3N4 
whose index of refraction is n = 2 [45]. 

The quasimonochromatic positive frequency compo- 
nent of the quantized field associated with these guided 
modes (g) at frequency wo takes the form 
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where v, is the group velocity, and 9 is the propagation 
constant of the guided modes. In the first Born approxi- 
mation the dispersive interaction of the guided field with 
Na atoms trapped near the surface of the waveguide at 
positions {r , Zn }, detuned far from resonance, is defined 
by the scattering equation [24], 

EM (r,t) = EO (r,t) 
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where @” is the atomic polarizability operator of the 
nt” atom, and 
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is the dyadic Green’s function for a dipole to radiate into 
the forward-propagating guided mode. In principle the 
Green’s function for an N,4-atom chain decomposes into 
a collective sub- and superradiant normal modes [23, 46], 
but in the far-detuning limit, these all are equally excited. 
The result is equivalent to the symmetric mode of inde- 
pendently radiating dipoles. The input-output relation 
for the mode operators then reads [24] 


FIG. 1. 
waveguides. (a) Electric field components of the H-polarized 


Fundamental guided modes of the nanophotonic 


HEi; mode of a circular nanofiber. From left to right are 
Re[uz(r_)], Re[uy(r_)] and Im[uz(r_)] in the ry-plane. (b) 
Same as (a) but for the H-polarized quasi-TEo; mode of a 
square waveguide. Black lines outline the waveguide bound- 
ary. The color scale is normalized to the maximum value 
of all field components for each waveguide mode. All other 
mode components not shown for both waveguide geometries 
vanish everywhere. (c) The normalized intensity distribution 
on the transverse plane for both geometries. The blue ar- 
rows indicate the local electric field’s direction and ampli- 
tude (relative length) at positions along the vertical waveg- 
uide axis, which only have an z-component of the mode. The 
stars indicate typical positions of trapped atoms (r‘, /a = 1.8 
for nanofiber [18] and similar scale for the square waveguide, 
r’, /w = 1.0, where a and w are the radius and width of the 
waveguides respectively). Dotted light gray lines show the 
corresponding V-mode contour which is the H-mode rotated 
by 90° around the waveguide propagation axis. The atom’s 
azimuthal position is chosen to be at a position with the V- 
mode being strongest. 


where 


Na 
bpp! = 2 en (5) 
is the phase operator associated with scattering polar- 
ization p' + p by a collective atomic operator. When 
p =p’, this is a phase shift; for p Æ p’, this leads to a 
transformation of the polarization of the guided mode. 
The Faraday effect arises from the irreducible rank- 
1 (vector) component of the polarizability tensor [47]. 
Given an atom with hyperfine spin f, this contribution 
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FIG. 2. (a) Schematic polarization spectroscopy geometry for 
the QND measurement and spin squeezing generation based 
on the Faraday effect. Atoms trapped near the surface of 
the nanophotonic waveguide cause a Faraday rotation of the 
guided light, which is measured in a polarimeter that detects 
the S2 component of the Stokes vector (intensity in the di- 
agonal D minus anti-diagonal D modes). (b) The evolution 
of the light’s polarization state on the Poincaré sphere (left 
to right). The Stokes vector of the light is prepared along 
the Sı direction, and the Faraday interaction causes a rota- 
tion around the S3-axis. Shot noise, shown as uncertainty 
bubble, limits the resolution of the detection. (c) The evolu- 
tion of the collective state before and after the measurement 
(left to right). The spin is prepared in a coherent state with 
projection noise shown as an uncertainty bubble. After the 
measurement the uncertainty in F, is squeezed, and the di- 
rection is correlated with the measurement outcome on the 
polarimeter. 


the D1- and D2-line transitions, respectively. We take the 
detuning, A, large compared to the excited state hyper- 
fine splitting. The resonant scattering cross section on a 
unit oscillator strength is co = 67/k2, where ko = wo/c. 
The polarization transformation associated with scatter- 
ing from H to V mode is determined by the operator 


dvi = i27 =a; [up (r1) x u(i) Ê, (6) 
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where F = r Ê) is the collective spin of the atomic 
ensemble. Thus, 
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and similarly for scattering from V to H. 
The polarization transformation can be expressed as a 
rotation of the Stokes vector of the light on the Poincaré 
sphere with operator components 
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By measuring the output Stokes vector in a polarimeter, 
we perform a QND measurement of a collective atomic 
operator to which it was entangled. In a proper con- 
figuration, this leads to squeezing of a collective spin. 


Launching H-polarized light corresponds to the initial 
Stokes vector along S1, and Faraday rotation leads to 
an S component, which is measured in a polarimeter 
(Fig. 2(a)). Taking the H-mode as a coherent state with 
amplitude p, the signal of the polarimeter measures 
Sout = (Byale™ + Bx,aeut)/2. Expressed in this way 
we see that the polarimeter acts as a homodyne detec- 
tor, with the input H-mode acting as the local oscillator 
and the photons scattered into the V-mode as the sig- 
nal. Formally, the input-output relation follows from the 
scattering equation, Eq. (7), and reads 


Sout = 53" +i(dvu—duv) SP =87 +xs(r)Ê Si. (9) 


The first term sin represents the shot-noise that funda- 
mentally limits the resolution of the measurement and 
thus the spin squeezing that can be obtained in a given 
time interval. The second term is the homodyne signal, 
where we have expressed the rotation angle around the 
3-axis of the Poincaré sphere as 
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We emphasize here the dependence of the rotation angle 
on the position of the atom in the transverse plane, r’, 
assumed equal for all atoms in the chain. In particular 
x3(r’,) depends on the overlap of uy(r’,) and uy(r’,) , 
indicating atomic scattering of photons from the H to V 
modes associated with the Faraday interaction. We have 
characterized this overlap by an effective area that defines 
the Faraday interaction at the position of the atom, 
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where ng = c/vg is the group index. A more tightly con- 
fined (smaller) area corresponds to a stronger interaction. 

By monitoring the Faraday rotation, we can perform 
a continuous measurement on the collective spin projec- 
tion F,. The “measurement strength,” which character- 
izes the rate at which we gain information and thereby 
squeeze the spin, is given by 
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where P;,, is the input power transported into the guided 
mode. The measurement strength is the rate at which 
photons are scattered from the guided H to V mode. De- 
coherence arises due to diffuse scattering into unguided 
modes and the accompanied optical pumping of the spin. 
In principle the photon scattering rate into 47 steradi- 
ans is modified over free space due to the Purcell effect, 
but we neglect this correction here. In the case of the 
nanofiber, this is a small effect at typical distances at 
which the atom is trapped [48, 49]. For the square waveg- 
uide, we will examine this correction in future work. 


Decoherence is due to optical pumping of the spin be- 
tween different magnetic sublevels. Henceforth, we re- 
strict to alkali atoms driven on the D1 or D2 line, with 
optical pumping rate 
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Here o (A) = 00 L is the photon scattering cross-section 
at the detuning A for a unit oscillator strength tran- 
sition, and the factor of 2/9 reflects the branching ra- 
tio for absorbing a m-polarized laser photon followed by 
spontaneous emission of a photon, causing optical pump- 
ing to another spin state. Im(r) = ngPnlug (ri)? = 
Pn/Ain(r’ ) is the input intensity into the guided H- 
mode at the position of the atom, where we have defined 
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to be the effective area associated with the input mode. 
We thus define the cooperativity /atom 
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This is our central result. Roughly, 1/[Arar(r’,)]?_ ~ 
Juy (r )|? luz (r )|?, thus C1 (r^) ~ ooluy(r’,)|?. In the 
context of homodyne measurement, the signal to be mea- 
sured is proportional to the overlap between the H- and 
V-modes, while the decoherence rate depends on the in- 
tensity of local oscillator or H-mode. How large the ini- 
tially unoccupied V-mode is at the atoms’ position deter- 
mines the signal-to-noise ratio for a QND measurement. 
We thus enhance the cooperativity by choosing the po- 
sition of the atoms so that the orthogonal, unoccupied 
mode is large, while the intensity of the local input mode 
that causes decoherence is small. 

We contrast this with squeezing arising from a birefrin- 
gence interaction, as we studied previously in [24]. Linear 
birefringence corresponds to a relative phase between the 
ordinary and extra-ordinary linear polarizations, which 
can arise both due to the geometry of the anisotropic 
modes relative to the placement of the atoms, as well as 
the atoms’ tensor polarizability. Here, the coupling is not 
optimal at the position of minimum intensity; it is max- 
imum at the angle 45° between the H and V modes. As 
such, one will not see as strong of an enhancement of the 
cooperativity as we will find in our protocol employing 
the Faraday effect. 

Figs. 3 and 4 show plots of 1/Apar, 1/Ain, and C1 
as a function of the position of the atom in the trans- 
verse plane, r'|, for the two nanophotonic geometries. 
We see that Apa, is essentially cylindrically symmetric 
sufficiently far from the surface for both the nanofiber 
and square waveguide geometries and thus the measure- 
ment strength is basically independent of the azimuthal 
position of the atoms. In contrast, Ain is azimuthally 
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FIG. 3. Contour plots of the effective mode areas and the 
cooperativity per atom near an optical nanofiber. (a) and (b) 
show the contour of the reciprocal effective Faraday inter- 
action mode area, Eq. (11), and the reciprocal input mode 
area in the xy-plane, Eq. (14), respectively. An x-polarized 
incident mode is assumed. (c) shows a contour plot of the co- 
operativity Eq. (15) in the xy-plane. The isovalue lines of Cy 
increase by 0.002428 on each gradient step from the outside 
inwards. The x and y coordinates are scaled in units of a for 
all three plots. 
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FIG. 4. Similar to Fig. 3 but for the square waveguide. In 
(a), the contour lines outside of the waveguide are essentially 
concentric circles. There are distortions near the four corners 
of the square waveguide shown in the plot mainly caused by 
numerical divergences. In (c), the isovalue lines of C4 increase 
by 0.005109 on each gradient step from the outside inwards. 
The x and y coordinates are scaled in units of w for all three 
plots. 


anisotropic. For input 2-polarization, 1/Ajn is the small- 
est along the y-axis at a given radial distance, which cor- 
responds to the smallest intensity of the input H-mode, 
and thus smallest optical pumping rate Yop. This angle 
corresponds to the position at which |uy(r’, )| is largest 
and thus yields the largest enhancement of C1. Thus, 
counterintuitively, we enhance the cooperativity by plac- 
ing the atom at the angle of minimum input intensity. 


This enhancement is even greater for the square waveg- 
uide which has more anisotropic guided modes compared 
to the cylindrical nanofiber. For typical geometries, given 
a nanofiber with radius a = 225nm and atoms trapped 
on the y—axis, a distance 200nm (0.8a) from the sur- 
face, the single atom cooperativity is Cı = 0.00728 at 
the optimal trapping angle; for the square waveguide of 
width w = 300nm, and atoms trapped 150nm from the 
surface, C1 = 0.0102 at optimum. Thus, with order 1000 
trapped atoms the N4-atom cooperativity is of order of 
10, sufficient to generate substantial spin squeezing. 


III. SPIN SQUEEZING DYNAMICS 


Given an ensemble of N4 atoms initially prepared in 
a spin coherent state for the hyperfine spin f, polarized 
in the transverse plane, e.g., along the x-axis, a QND 
measurement of the collective spin F, will squeeze the 
uncertainty of that component. The metrologically rele- 
vant squeezing parameter is [50], 


p = Nal(AF?) 
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Under the assumption that the state is symmetric with 
respect to exchange of any two atoms, valid when we start 
in a spin coherent state and all couplings are uniform 
over the ensemble, the collective expectation value can 
be decomposed into 
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The first term of Eq. (17a) and Eq. (17b) is the projec- 
tion noise associated with the N4 identical spin- f atoms, 
and the second term of Eq. (17a) is determined by two- 


body covariances, (APAP, = (APPA) = 
ifj 
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(i?) = (FEN Fy, Negative values in these two- 
body correlations correspond to the pairwise entangle- 
ment between atoms, leading to spin squeezing [51]. 
Note, when the detuning is sufficiently far-off resonance, 
all collective sub- and super-radiant modes [23, 46] are 
equally (and thus symmetrically) excited. In this paper, 
we work in the dispersive regime with a few thousands of 
atoms and can safely ignore the atom-atom interaction 
caused by multiple scattering, and hence the collective 
atomic system satisfy the exchange symmetry. 

To study the spin squeezing dynamics, we follow the 
method first developed by Norris [52]. We employ a first- 
principles stochastic master equation for the collective 
state of N4 atoms, 


dô = dôļlonp T dôlop $ (18) 


The first term on the right-hand side of Eq.(18) governs 
the spin dynamics arising from QND measurement [9, 


53], 
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where « is the measurement strength defined in Eq. (12), 
and dW is a stochastic Wiener interval. The conditional 
dynamics are generated by superoperators that depend 
on the collective spin 


H [ô] = Fp + pF. — UF.) p, 
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The second term governs decoherence arising from op- 
tical pumping, which acts locally on each atom, dplop = 
EA D(”) [A] dt, where 
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Here Ñ (o) is the effective nonHermitian Hamiltonian de- 
scribing the local light shift and absorption by the nt” 


atom and Wo” is the jump operator corresponding to 
optical pumping through absorption of a laser photon 
followed by spontaneous emission of a photon of polar- 
ization q [47] (see Appendix A). 

The rate of decoherence is characterized by the optical 
pumping rate, Yop. Note, the optical pumping superoper- 
ator, Eq.(21), is not trace preserving when restricted to a 
given ground-state hyperfine manifold f. Optical pump- 
ing that transfers atoms to the other hyperfine mani- 
fold in the ground-electronic state is thus treated as loss. 
Moreover, if the atoms are placed at the optimal posi- 
tion in the transverse plane, the local field of the guided 
mode is linearly polarized. In that case the vector light 
shift vanishes, and for detunings large compared to the 
excited-state hyperfine splitting, the rank-2 tensor light 
shift is negligible. The light shift is thus dominated by 
the scalar component, which has no effect on the spin 
dynamics. In that case Aug = — (Meri, representing an 
equal rate of absorption for all magnetic sublevels. 

Following the work of Norris [52], the solution to the 
master equation is made possible by three approxima- 
tions. Firstly, we restrict the subspace of internal mag- 
netic sublevels that participate in the dynamics. The 
system is initialized in a spin coherent state, with all 
atoms spin-polarized along the x-axis. We denote this 
as the “fiducial state” | = |f, Mme = f). Through QND 
measurement, spin squeezing is induced by entanglement 
with the “coupled state” |4} = |f, Me = f-—1). Op- 
tical pumping is dominated by “spin flips” M —> |J) 
and “loss” due to pumping to the other hyperfine level. 
We finally include a third internal magnetic sublevel 
|T} = |f, mz = f — 2) to account for “transfer of coher- 
ences” that can occur in spontaneous emission [52, 54] 
(see Fig. 5). Restricted to this qutrit basis, the internal 


hyperfine spin operators are 
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where we have defined the atomic population and coher- 
ence operators Gyq = |b) (al. 
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FIG. 5. Schematic energy level diagram for cesium atoms 
probed on the D2 line, 6S1/2 — 6P3/2. Relevant dynam- 
ics are restricted to a truncated qutrit subspace of ground 
levels. Atoms are prepared in the fiducial state |t), and 
driven by x-polarized (7) light. The Faraday rotation cor- 
responds to coherent scattering of 7 — o in this basis, and 
measurement backaction leads to entanglement between pairs 
of atoms in |t) = |f =4,mz = 4) and the coupled state 
|L) = |f =4,mz = 3). Optical pumping can cause spin flips 
It) > H) and H} > |T} = |f =4,mz = 2). The latter pro- 
cess is included to account for transfer of coherences. The 
detuning A is taken to be large compared to the excited state 
hyperfine splitting. 


Secondly, we assume the collective state is symmetric 
under exchange of spins. This approximation is valid 
when all atoms see the same probe intensity and in 
the far-detuning regime when all sub- and super-radiant 
modes are equally excited. With this, we can limit our 
attention to the symmetric subspace and define, for ex- 
ample, the symmetric two-body covariances by 
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[Ao Aa’) ARA), (23) 
where the superscripts, (1) and (2), label arbitrary 
two different atoms in the ensemble. Due to the ex- 
change symmetry, (Ag Ao). = (Ad AP) = 
(Ac? AoW Js which reduces the number of n-body mo- 
ments required to simulate the spin dynamics of the en- 
semble. 

Thirdly, we make the Gaussian approximation, valid 
for large atomic ensembles, so that the many-body state 
is fully characterized by one- and two-body correlations. 


Equivalently, the state is defined by the one- and two- 
(1) 


body density operators, with matrix elements pġ = 


(Ge); pea, = (Ac Ao?) | in the symmetric sub- 
space. We track the evolution of the correlation func- 
tions through a set of coupled differential equations [52]. 
Optical pumping, acting locally, couples only n-body cor- 
relations to themselves, e.g., 


d(Qoyq Mog. ) 


Slop 


= (Di [Aol Ao?) dt + (AoW) Dt [Ao )) dt. (24) 
QND measurement generates higher order correlations 
according to 


. K 
d Sea) lonp ~ | 


(Lİ [Sta] yar (HÌ [Soa] )\dW. (25) 
We can truncate this hierarchy in the Gaussian approxi- 
mation, setting third order cumulants to zero. Thus, for 
example, 

d(Aoyg Arge ) 


S 


QND 
QND (50a) oND (4a) lan) 


= Caolonp (4@u2)| anv) z dora) onp d(oac) ene 
= —K(AofP AF.) (AF,Ao®) dt, (26) 


= dha On), 


where we have employed the Ito calculus dW? = dt. 
Note, by setting the third-order cumulants to zero, the 
contribution of the £ superoperator to dynamics of 


the two-body covariances vanishes, d(Aol? Ao®)) 


Sle 
s(cÌ [AR Ag] ),dt = 0. This indicates the events 
of no-photon detection under the Gaussian state approx- 
imation do not affect atom-atom correlations; the mea- 
surement backaction and squeezing arise from the homo- 
dyne detection in the guided modes. 

Using all of the approximations above, the collective 
spin dynamics can be efficiently calculated for the ensem- 
ble of qutrits (dimension d = 3) with d? = 9 equations 
for the one-body quantity, (va), and d?(d? + 1)/2 = 45 
equations for the two-body covariances, (AG Kal?’ 
in the symmetric subspace independent of the number of 
atoms. With this formalism in hand, we can calculate 
the squeezing parameter, Eq.(16), as a function of time 
by finding time-dependent solutions for the one-body av- 
erages ( fa) and (A i), and the two-body covariances 


(AFP AFP). The detailed approach to calculating the 
collective spin dynamics can be found in Appendix A. 
Using this formalism, we calculate the squeezing of an 
ensemble of cesium atoms, initially spin-polarized in the 
6S; /2,|f = 4, Mg = 4) state. We choose the guided mode 
frequency near the D2 resonance, 661/2 — 6P3/2, de- 
tuned far compared to the excited-state hyperfine split- 
ting. In Fig. 6, we plot the reciprocal of the spin squeez- 
ing parameter as a function of time and its peak as a func- 
tion of atom number, N4, for both the optical nanofiber 


and the square waveguide cases. By placing atoms 200 
nm away from the nanofiber surface (r’) = 1.8a), our 
simulations for 2500 atoms yield 6.3dB of squeezing. Us- 
ing the square waveguide platform with the same num- 
ber of atoms placed 150nm away from the surface, our 
calculation yields 12.9dB squeezing. As we have shown 
in Sec. II, the square waveguide geometry enhances the 
anisotropic contrast of the two orthogonal guided modes 
and dramatically reduces the relative local intensity when 
the atoms are placed on the y-axis. This results in a 
large cooperativity, higher peak spin squeezing, achieved 
in shorter time, and with a relatively slower decay when 
compared to the nanofiber. In addition, in Figs. 6(b) 
and 6(d) we show how the peak squeezing scales with 
the number of trapped atoms when the atom positions 
are fixed as above. 
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FIG. 6. The reciprocal spin-squeezing parameter, Eq. (16). 
Subfigs. (a) and (c) are plots of €~? as a function of time 
in units of the optical pumping rate yop, for the cylindrical 
nanofiber and square waveguide, respectively, for Na = 2500, 
with other parameters given in the text. These curves peak 
at time determined by the detailed balance of reduced un- 
certainty due to QND measurement and decoherence due to 
optical pumping. Subfigs. (b) and (d) show the plots of the 
peak value €~? as a function of N4 for the nanofiber and 
square waveguide correspondingly. 


In the absence of any other noise, the cooperativity 
of atom-light coupling increases as the atoms are placed 
closer to the waveguide surface. Figs.7(a) and 7(c) show 
the peak squeezing as a function of rı for both nanofiber 
and square waveguide geometries with 2500 atoms. With 
the same setting, we also plot out the cooperativity, C1, 
on a logarithm scale in Figs.7(b) and 7(d) as a function 
of atom radial distance to the center of both waveguide 
geometries. We find that C, is proportional to e~°" and 
the peak squeezing scales as VOD, where 6 ~ 1.65/a 


for the nanofiber and 6 ~ 2.14 x 2/w for the square 
waveguide with 2500 atoms. The cooperativity of the 
square waveguide increases faster than the nanofiber as 
the atoms approach to the waveguide surface. 
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FIG. 7. Peak squeezing parameter ((a) and (c)) and coopera- 
tivity at the optimal azimuthal trapping position ((b) and (d)) 
as a function of the radial distance to the waveguide axis, for 
Na = 2500, with other parameters given in the text. (a) 
and (b) are for the nanofiber case; (c) and (d) are for the 
square waveguide case. 


IV. SUMMARY AND OUTLOOK 


In this paper we studied the cooperativity of the atom- 
photon interface for two nanophotonic geometries: a 
cylindrical nanofiber and a square waveguide. Due to the 
anisotropic nature of the guided modes, one can strongly 
enhance the cooperativity by trapping atoms at positions 
that maximize the rate at which photons are forward- 
scattered into the orthogonally polarized guided mode, 
while simultaneously minimizing the rate of which they 
are scattered into free space. Counterintuitively, the op- 
timal geometry is such that atoms at a certain distance 
from the surface are trapped at the azimuthal angle of 
minimal intensity of the probe. We applied this idea to 
study the generation of a spin squeezed state of an ensem- 
ble atoms induced by QND measurement, mediated by 
the Faraday interaction in the dispersive regime. With 
realistic parameters, our simulation shows more than 6dB 
of squeezing for a cylindrical nanofiber or 12dB for the 
square waveguide with 2500 atoms. The amount of spin 
squeezing we predict based on Faraday effect is substan- 
tially larger than for the birefringence-based spin squeez- 
ing protocol studied in our earlier work [24]. Although 
we have only considered a cylindrical nanofiber and a 


square waveguide geometries, the ideas presented in this 
paper are applicable to other nanophotonic waveguide 
geometries which could show enhanced cooperativity. In 
addition, our model of decoherence, simplified when the 
detuning is large compared to the excited state hyperfine 
splitting, is almost certainly not the optimal operating 
condition. 


Our simulations are based on a first-principles stochas- 
tic master equation that includes QND measurement- 
backaction that generates entanglement and results in 
spin squeezing, as well as decoherence due to optical 
pumping by spontaneous photon scattering [9, 24, 52]. 
This simulation is made possible by a set of simplifying 
assumptions: (i) we restrict each atom to a qutrit, em- 
bedded in the hyperfine manifold of magnetic sublevels; 
(ii) the state is exchange-symmetric with respect to any 
two atomic spins; (iii) the many-body state is fully char- 
acterized by one- and two-body correlations (the Gaus- 
sian approximation). With these, we solve for the metro- 
logically relevant squeezing parameter as a function of 
time and see the tradeoffs between QND measurement 
and decoherence for various geometries and choices of pa- 
rameters. Our method is extendable to include higher or- 
der correlations, which become manifest for large squeez- 
ing, when the Holstein-Primakoff (Gaussian) approxima- 
tion breaks down. The computational framework we have 
developed here should allow us to study n-body moments 
with an acceptable computational load. 


In future work we intend to extend our analysis in 
a number of directions. While we focussed here on 
the enhancement of the cooperativity in a nanophotonic 
waveguide-based QND measurement, we did not fully an- 
alyze the impact of Purcell effect and the modification of 
spontaneous emission rates in our simulations. We have 
also neglected here the motion of atoms in the optical 
lattice. In practice, however, these effects are impor- 
tant, the latter having been observed in the nanofiber 
experiments [37, 38, 55, 56]. We will include these in 
future studies, with an eye towards the development of 
new strategies for atomic cooling and state initialization 
in the nanophotonic platforms [57]. We expect that our 
proposed geometry, which places the atoms at positions 
of minimum intensity, could also help reduce the pertur- 
bation on the motion of trapped atoms due to the probe 
that causes a disturbance in the signal [37]. 


Finally, we have studied here the dispersive regime of a 
QND measurement, where the probe light is detuned far- 
off-resonance, and multiple scattering of photons among 
atoms is negligible. To extend our theory, it is necessary 
to include the collective effects such as super and sub- 
radiance [23, 46, 58], with applications including quan- 
tum memories [34, 59], the generation of matrix product 
states and cluster states [60—64]. 
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Appendix A: Modeling collective spin dynamics 


In this Appendix we give further details of the equa- 
tions of motion for spin squeezing as a function of time, 
as discussed in the main text, building on the work of 
Norris [52]. In the symmetric subspace, and in the Gaus- 
sian approximation, we track a set of one- and two-body 
correlation functions that determine the metrologically 
relevant squeezing parameter, Eqs. (16,17). The evolu- 
tion is determined by the dynamics induced by the QND 
measurement and decoherence due to optical pumping, 
Eqs. (18-21). For concreteness, we consider here ce- 
sium atoms, initially prepared in a spin coherent state, 
with all atoms in the stretched state polarized along 
the x-axis, |t) = |651/2, f = 4, Mı = 4). The atoms are 
trapped on the y-axis at a distance r’, from the core 
axis of the waveguide and are probed with guided light 
in the H-mode, which has linear polarization in the x- 
direction at the position of the atoms. We take the 
detuning large compared to the excited state hyperfine 
splitting, e.g., 4 GHz red detuned from the D2 line, 
1651/2, f =4) > |6P3/2, f’ = 3). Spontaneous emission 
from the probe may result in optically pumping of atoms 
to the other hyperfine manifold, f = 3; we treat this as 
a loss channel under the approximation that, over the 
time interval of interest, there is a negligible probability 
for these atoms to repump to f = 4. We also include a 
bias static magnetic field along the z-axis. This does not 
affect the QND measurement of F,, but ultimately, we 
must calculate all dynamics in the rotating frame. 

Spontaneous emission, optical pumping, and the re- 
sulting decoherence acts locally on each atomic spin, gov- 
erned by the master equation, Eq. (21). For light linearly 
polarized in the x-direction, and for detunings large com- 
pared to the excited state hyperfine splitting, this takes 
the simplified form [47] 


dp 


= Dio” 
dt [op 

es AG a(n), Yoo fmen, FMaM) PO 
=- lA Aron + Ta (FP OOF HAL OFP), (AL) 
where FO = Ain) ti FO are the raising and lowering 
operators for projection of spin along the z-axis, and 

= ees lin + : : . 
Yoo = TAF Top is the optical pumping for linear po- 


larization in the far-detuned limit, given intensity Jj, 
at the position of the atom (we assume that all atoms 
are trapped the same distance for the waveguide, and 


thus see the same intensity). The atomic Hamiltonian is 


Hy, = vn AQ fs” + Hyg, the sum of the Zeeman inter- 
action with a bias magnetic field, giving rise to Larmor 
precession at frequency Qo, and the light shift, Hrs = 
hx3F.53, due to the probe. Even for far detuning, the 
residual tensor light shift cannot be neglected as it scales 
at 1/A?, the same as Yop and «. In principle a two-color 
probe can remove the tensor term which would otherwise 
lead to a degradation of the mean spin, and thus a reduc- 
tion metrologically useful squeezing [65, 66]. We neglect 
this effect here. Finally, going to the rotating frame at the 
Larmor frequency, f{”) — f cos(Qot) + FL” sin(Qot), 
fi > A cos(Qot) — ae sin(Qot), and averaging the 
rapidly oscillating terms (RWA), the master equation 
takes the form 


dp 
dt 


=> —YopP 
an (Aep OE ARa Pp i, (A2) 


With atoms initially prepared in the “fiducial state,” 
It) = |f =4,m, = 4), we include in the dynamics the 
“coupled state,” |J) = |f = 4, Mmg = 3), and the “trans- 
fer state,” |T) = |f =4,m, = 2). Restricted to this 
qutrit substate, the spin projector operators { A Fis I 
are given in Eqs. (22). 

With all the components of the stochastic master equa- 
tion defined in Eq. (18-21), one can derive the equa- 
tions of motion of one- and two-body moments straight- 
forwardly using the symmetric Gaussian state approxi- 
mation. Some explicit results have been given in Eq. (24 
-26). The equation of motion for the optical pumping 
dynamics of the one-body moment (Sba), is given by 


(6a) lop = (D [Goa] dt =) tr (D! [Gba]ac)Gac)dt. (A3) 
c,d 


The equations of two-body moments of the optical pump- 


ing can be derived similarly. Continued from Eq. (24), 
we have 


daot A02) | 
op 


(ADİ [a] Ag) dt + (Aoi? ADI GA )s dt 
= S-tr( Di [őba]őmn) (Ao®, A02), 


min 


ps 


In deriving Eqs. (25) and (26), we have used the Gaus- 
sian state assumption to write three-body moments in 
connection to one- and two-body moments, (ABC) = 


(4B)(C) + (AC)(B) + (BC)(A) - 2(4)(B)(C). 1 


[Fae ]émn) (Aee D Aon), (A4) 
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only keeping coherence operators of the nearest coupling 
states, Eqs. (A3) and (A4) recover the same results found 
by Norris [52]. 

We apply similar techniques to derive the equations of 
motion due to QND measurement Eqs. (25) and (26) to 
yield 


dG ba) au = 


2 tr (LT [Goa] Fac) (Fac) dt 
c,d 


k DI. (H' [Spaleac) GacydW. (A5) 


c,d 


QND 


=; v (at (601) +551 Gra) + 5a) Gor) +5or(F1a)) 
: T= (5asGor)+ or Gra) t Sar Goy)+ sural 


2f-1 


-VE er Crt en- A ud (Our) +r) 
toran [yf (aot? aoao) 


[of—1 
+ a (Aog Aot, J + (Aona Ao l 5) | 


i (as) +a, Öre) +5e1 Car) + at Öre) 


f 


j > Oe, Gar)+ bar Oe) Her Gay) + 5ay 6 n) 


- Jhe ente- oad enter) 


+(N4-1) EC Act) Acs?) + (Act? Ao{?),) 


H= al = (oft Ao) (Ach Aog.) 


By combining the optical pumping and QND measure- 
ment contribution to Eqs. (A3-A6), one can find a set of 
stochastic equations of a closed set of variables of 


{( Gra); (Aba Aca), la, b, c,d € MTh} (A7) 


Jar (A6) 


in the symmetric qutrit subspace. The matrix of equa- 
tions is sparse and close to diagonal which indicates 
only nearest neighboring coupling is possible in the 
{ (va), (Aba Aac), } basis. In the symmetric qutrit 
subspace, we have 45 two-body moment variables and 
corresponding sparse second-order equations, which we 
solve numerically. 


